Adaptive  Observers  for  Structurally  Perturbed  Positive  Real  Delay  Systems 


M.  A.  Demetriou*  R.  F.  Curtain  t  K.  Ito* 


ABSTRACT 

The  aim  of  this  investigation  is  to  construct  an  adaptive  observer  for  a  class 
of  multivariable  delay  systems  having  a  known  exogenous  input  and  a  struc¬ 
tured  perturbation  with  an  unknown  constant  parameter,  such  as  the  case  of 
static  output  feedback  with  an  unknown  gain.  The  adaptive  observer  uses  the 
nominal  dynamics  of  the  unperturbed  plant  and  an  adaptation  law  based  on 
the  Lyapunov  redesign  method.  We  obtain  conditions  on  the  system  to  ensure 
uniform  boundedness  of  the  estimator  dynamics  and  the  parameter  estimates, 
and  convergence  of  the  estimator  error.  We  illustrate  our  approach  with  a 
multivariable  delay  example  for  which  we  obtain  convincing  numerical  results. 


1  Introduction 


In  Curtain,  Demetriou  and  Ito  [2]  we  constructed  adaptive  observers  for  the  following  class 
of  infinite-dimensional  systems  on  a  given  Hilbert  space  X  with  inner  product  and  norm 
{•,  •)  and  |  •  | Xi  respectively: 

=  Apx(t)  +  Bu(t)  +  /(<);  x(0)  =  xq  €  X  (1.1) 

y{t)  =  Cx(t),  (1.2) 

where 

AP  =  A0  +  BTC  (1.3) 

and  Ao  is  a  generator  of  an  exponentially  stable  Co-semigroup  T(/,),  i.e.  there  exists 
constants  M,  /i  >  0  such  that 

\\T(t)\\  <  Me^‘.  (1.4) 

The  signals  u(t)  and  y(t)  are  the  vector-valued  inputs  and  outputs,  respectively.  f(t)  is 
an  X- valued  known  exogenous  input  and  B  €  C ( (D  .X).  C  €  C{X.  (D ) .  Consequently, 
Ap  generates  a  strongly  continuous  Co-semigroup  for  any  T  and  the  system  is  well-posed. 
A  key  assumption  is  that  the  original  system  (Ao,B,C)  satisfy  a  positive-real  condition, 
where  by  positive-real  we  mean  the  following. 


Definition  1.1  Suppose  that  G(-)  :  (Uq  — >  C((Um)  where  (Pq  =  {s  :  Res  >0}.  If 


1A- 

ef 

ii 

03 

(1.5) 

G(s)  is  holomorphic  on  (Dq 

(1.6) 

G(s)  +  G(s)*  >  0  for  all  s  =  jcu ,  lo  E  i?, 

(1.7) 

then  G  is  positive- real. 
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In  fact,  we  require  that  G(s  —  y)  —  C  ((s  —  y)I  —  Aq) B  be  positive-real  for  some  y  >  0. 

The  motivation  is  that  the  perturbation  BTC  arises  as  the  unknown  constant  gain  in 
the  feedback  law  u  =  Ty.  Thus,  the  actual  formulation  is 

x(t)  —  A0x(t)  +  BTy(t)  +  Bu(t)  +  /(<)*  x(0)  =  x0-  (1.8) 

The  proposed  state  estimator  (observer)  is 

x(t)  =  A0x(t)  +  Bu(t)  +  BT(t)y(t)  +  f(t),  2(0)  =  x0,  (1.9) 

where  x(t)  is  the  state  estimate  at  time  t  and  T ( /; )  is  the  adaptive  estimate  of  the  unknown 
gain.  The  adaptation  rule  for  T(f)  was  based  on  the  Lyapunov  redesign  method  [5,  6]  that 
has  proved  successful  for  finite-dimensional  systems  and  it  depends  on  the  existence  of 
a  solution  to  a  certain  constrained  Lyapunov  equation.  The  following  is  the  main  result 
from  Curtain,  Demetriou  and  Ito  [2], 


Theorem  1.2  Consider  the  structurally  perturbed  system  (l.l)-(l.S),  where  A$  is  the 
generator  of  an  exponentially  stable  semigroup  on  X,  B  E  C(Cm ,  X) ,  C  E  C(X.  (Um),  f(t) 
is  a  known  exogenous  signal  with  values  in  X,  and  T  is  an  unknown  matrix  feedback  gain. 
If  there  exist  Q  E  C(X)  and  L  E  C(D(A),X)  or  £(D(A),  (Tm)  satisfying  the  constrained 
Lyapunov  equation  for  x  E  D(A) 

(Ao  +  yl)*  Qx  +  Q  {Aq  +  yl)  x  =  —L*Lx  (1.10) 

B*Qx  =  Cx ,  (1.11) 


then  the  state  estimator  defined  by  (1.9)  and  the  adaptation  ride  with  adaptation  matrix 
gain  G  =  GT  >  0  given  by 


f  (<)  =  GCe(t)yT(t) 

f  (0)  =  f0 


(1.12) 


have  the  following  properties. 


(i)  If  u,y  E  1/00(0,00;  (Um),  then  r(f),  x(t)  and  the  estimation  error  eft)  =  x(t)  —  x(t) 

are  bounded  in  norm  for  t  >  0  and  \\e^tQ^  (x(t)  —  x(t))  \\x  — ¥  0  as  t  — >  oo. 

(ii)  I  If  y  E  Loo(0,  oo;  C  m )  fl  1/2(0,  oo;  C  m)  then  ||e 2"*  (x(t)  —  x(t))  ||x  — t  0  as  t  — >■  oo 


(in)  If  the  conditions  in  (ii)  hold  and,  in  addition,  the  plant  is  persistently  exciting,  i.e. 
there  exists  To,  5o  and  cq  such  that  for  each  sufficiently  large  t,  >  0  there  exists 
t  E  [t.t  +  Tq]  such  that 


>  eo 


for  every  unit  vector  w  E  Rm ,  then  we  can  achieve  parameter  convergence,  i.e. 


T(t)  — >■  T,  as  t  — >•  oo. 


The  crucial  assumption,  that  the  constrained  Lyapunov  equation  (1.10)-  (1.11)  has  a 
solution,  is  a  difficult  problem  for  infinite-dimensional  systems  and  only  partial  results  are 
available.  In  Curtain,  Demetriou  and  Ito  [2],  we  designed  adaptive  observers  for  scalar 
systems  of  parabolic  and  delay  type  by  exploiting  recent  results  on  constrained  Lyapunov 
equations  in  Curtain  [1],  While  the  theory  also  applied  to  multivariable  systems,  the 
conditions  prove  to  be  difficult  to  verify  in  the  examples.  It  is  the  aim  of  this  paper 
to  apply  the  theory  to  design  an  adaptive  observer  for  a  class  of  multivariable  retarded 
systems.  We  demonstrate  the  power  of  the  results  by  numerical  simulations. 


2  A  class  of  positive-real  retarded  systems 

We  consider  the  following  class  of  retarded  systems 


x(t)  =  Aox(t)  +  ax(t  —  h)  +  Bou(t ); 
a;(0)  =  r,  x(0)  =  f(0),  —h  <  6  <  0; 
y{t)  =  C0x(t), 


(2.1) 

(2.2) 

(2.3) 


where  h  >  0  represents  the  point  delay,  x(t)  €  (D  n,  Ao  €  £((D  n),Bo  €  £(<D  m,(D  "),Co  € 
£((D  ",<E  m),r  e  (D  n,  and  /  €  L2(— h,  0;  C  n).  The  delay  system  (2.1)-  (2.3)  can  be  for¬ 
mulated  on  the  state  space  X  —  (D  "  ®  L2(— h,  0;  (D  ")  with  generating  operators  defined 
by 

X  B0u  \  (  r  X 

o  ’  C  /(•) 


Bu  = 


=  Cor, 


(2.4) 


f  \ 

^  I  e  X|  /  is  absolutely  continuous, 
^(•)eL2(-fi,0;(D")and/(0)  =  r 


r 

/(•) 


A0r  +  af(-h)  \ 

•  (2-5) 

d6»  / 

A  generates  a  Co  semigroup  T(-)  on  X  and  B.  C  are  bounded  operators  (see  Curtain  and 
Zwart  [3],  chapter  2.4  and  Theorems  4.2.6,  4.2.10).  The  spectrum  of  A,cr(A)  consists  of 
eigenvalues  which  are  the  solutions  of  det  (A(A))  =  0,  where  for  A  €  (D 


A(A)  =\I-A0-  ae-Xh. 

T(-)  is  exponentially  stable  if 

cr(A)  C  {s|Re  s  <  — /i}  for  some  /i  >  0. 

The  transfer  function  G  of  the  system  is  given  by 

G(s)  =  CqA^Bq. 

We  also  need  to  consider  a  more  general  type  of  observation  operator 


Ce 


r 

/(•) 


=  Cir  +  C2f(-h), 


(2.6) 

(2.7) 

(2.8) 

(2.9) 


where  Cj  €  C ( C  n,(D  q).]  =  1,2.  The  transfer  function  H  of  the  system  defined  by  (2.1), 
(2.2)  and  the  observation 

y(t)  =  C\x(t)  +  C2x{t  -  h )  (2-10) 

is  given  by 

E{s)  =  {C1  +  C2e-hs)A{s)-1B0.  (2.11) 

We  now  give  conditions  under  which  the  system  (2.1)-(2.3)  will  be  positive-real. 

Lemma  2.1  Suppose  that  (Ao,Bo)  is  controllable,  Ao  +  \a\I  is  a  stable  semisimple  matrix 
and  Go(s)  =  Cq(sI  —  Ao)~1Bo  is  positive-real.  Then  there  exist  matrices  Lo,  and  P  = 
P*  >  0  such  that 


A*0P  +  PAo  =  -L*0L0 
B*0P  =  Co 

Moreover,  G(s)  given  by  (2.8)  is  positive-real. 


(2.12) 

(2.13) 


Proof.  (2.12),  (2.13)  follow  from  the  finite-dimensional  positive-real  results.  We  rewrite 
(2.12)  as  follows 


(A0  +  \a\I)*P  +  P{A0  +  \a\I)  =  -L*0L0  +  2|a|P  (2.14) 

Denoting  K  =  LqLq  —  2|a|P  and  A\  =  Aq  +  \a\I,  we  deduce  the  following  from  (2.14) 

jt{eMtPx,eMtPx)  =  -{KeAltPx,eAltPx ) 
and  since  A\  is  a  stable  matrix 

POO 

(Px,  x)  —  /  (KeAlt  Px,eAlt  Px)  dt. 

Jo 

Now  Ai  is  semisimple  and  so 


eAltx  =  ^ ~2eXit(x,ei)ei , 
i= l 


where  (Aj,  e()  are  the  eigenvalue-eigenvector  pairs  of  A\.  Thus 


(Peuei)  =  [ 

Jo 


2ReXit(Kehei)dt 
{ Kehei ) 


2 Re  A. 

and  Re  \  <  0  since  A\  is  stable,  and  hence 

K  =  L*0L0  -  2\a\P  >  0. 


An  easy  calculation  gives 
n(iw) 

=  G(iu)*  +  G(iw) 

=  B£{-iu  -  AI  -  ae^b)-1  (LqLq  -  2aP cos Luh)(iiu  -  A0  -  ae^by1  B0 

>  By-ico  -  Aq-  ae^by1  (LqLq  -  2\a\P)(ico  -  A0  -  ae^^Po 

>0, 

where  we  have  used  (2.12)-(2.14).  Finally,  A  is  stable  if  Ao  +  \a\I  is.  □ 

Next  we  obtain  a  spectral  factorization  for  II  of  the  form 

II(ia;)  =  H(uj)*H(«d).  (2-15) 

Lemma  2.2  There  exist  square  matrices  C\,C2  €  C((Tn)  such  that  (2.15)  holds  with  the 
spectral  factor  H  given  by  (2.11).  The  matrices  are  any  solutions  of 

CfC)  +  C)C2  =  C(C0:  C(C2  =  CqCh  =  —2aP.  (2.16) 

Proof.  If  there  exist  (71,(72  satisfying  (2.16)  it  is  straightforward  to  verify  that  S 
satisfies  (2.15).  So  it  remains  to  show  the  existence  of  Ci,C2.  Without  loss  of  generality 
we  can  take  C\ ,  C2  to  be  symmetric  and  a  to  be  positive.  We  choose  a  square  matrix  F 
such  that 

(Ci  +  C2)*(Ci  +  C2)  =  <7o(70  -  2 aP  =  F*F. 

Then,  substituting  C2  =  F  —  C\  we  obtain  the  Riccati  equation  for  C\: 


C\F  +  F*C\  -  2C\C\  +  4aP  =  0. 


This  has  a  solution  which  we  take  to  be  C\  and  Co  =  F  —  C\.  □ 

While  in  finite  dimensions  the  existence  of  a  spectral  factor  and  controllability  is  suf¬ 
ficient  to  ensure  the  existence  of  a  solution  to  the  constrained  Lyapunov  equation,  the 
infinite-dimensional  case  is  more  subtle.  We  refer  to  Curtain  [1],  section  7,  for  the  rest 
of  the  proof  of  the  existence  to  a  solution  to  (2.12),  (2.13).  Essentially,  it  amounts  to 
showing  that  the  observation  (2.9)  is  well-posed  and  this  is  known  to  be  the  case.  So  we 
may  conclude  that  Theorem  1.2  holds  under  the  assumptions 

•  (Ao^Bo)  is  controllable 

•  Ao  +  |  a  1 1  is  a  stable  matrix 

•  Gq(s)  =  Co(sI  —  Aq)~1Bo  is  positive-real. 


3  Examples  and  Numerical  Results 


We  used  the  following  matrices  in  (2.1),  (2.3)  to  simulate  the  proposed  adaptive  scheme 


Ao 


-1.0  0.4 

0.2  -2.0 


Bo  = 


1.00  0.01 
0.01  1.00 


Co 


2.01  1.04 
1.02  4.01 


The  pair  ( Ao-,Bo )  is  indeed  controllable  and  the  solution  to  the  constrained  Lyapunov 
equation  is 


P  = 


2  1 
1  4 


with  LqLq 


3.6  1.4 

1.4  15.2 


One  can  easily  check  that  BqP  =  C  and  thus  Co(sI  —  Aq )  1 B  is  positive  real.  The  scalar 
a  is  chosen  as  a  =  0.01  with  a  delay  h  =  1.  In  this  case  we  have  that  eig(To  +  |ck| I)  = 
(-0.8255,-1.9745)  which  gives  a  stable  Ao  +  \a\I  matrix.  Therefore  all  the  conditions 
required  for  Theorem  1.2  to  hold  are  satisfied. 

The  unknown  parameter  in  the  perturbation  term  is  taken  to  be 


'  Tn 

ri2 ' 

’  4.0  x  10~4 

1.0  x  10~4  ’ 

r21 

r*22 

2.0  x  10~4 

2.0  x  IQ-1 

The  initial  conditions  for  the  state  x.  its  estimate  x  and  the  parameter  T(f)  are  given  by 


x{9) 


sin(4#  —  1)  +  sin(l) 
cos  (46  —  1) 


x(6) 


0.6  (sin(4#  —  1)  +  sin(l)) 
1.4  (cos(40  —  1)) 


where  6  6  [—1,0]  and 


f(0) 


2.0  x  lO^1  0.5  x  10~4 
1.5  x  10~4  3.0  x  lO^1 


The  plant  input  for  this  set  of  simulations  was  taken  to  be 


u{t) 


cos(5.07t  t)  +  sin(7.07rt)  +  sin(24.27rf) 
sin(5.07rf)  +  sin(8.07rf)  +  cos(29.2-7rf) 


The  adaptive  laws  for  the  entries  of  the  gain  matrix  become 

?n(f)  =  -5  x  10~2ei(f)yi(f),  fi 2{t)  =  -3  x  10~5£i(t)y2{t), 
?2iW  =  -5  x  10~5£2(%i(t),  f22(t)  =  -5  x  10~2£2(%2(f), 


where  e^t)  =  y^t)  -  &(<),  i  =  1, 2. 


x1  (t)  (dashed)  and  its  estimate  (solid) 


Figure  1:  Evolution  of  (a)  x\  (dashed)  and  xi  (solid),  (b)  X2  (dashed)  and  X2  (solid). 


The  above  system  was  simulated  using  the  approximation  scheme  developed  in  [4]  with 
N  =  12  splines  to  discretize  the  interval  [—1,0],  thus  resulting  in  26  equations  for  the  state, 
26  for  its  estimate  and  4  for  the  parameter  matrix  T,  totaling  56  equations.  Due  to  the 
infinite  dimensionality  of  the  delay  system,  the  discretization  index  N  =  12  was  chosen  as 
the  smallest  possible  one  that  would  provide  a  numerically  reliable  approximation  to  the 
above  system.  The  resulting  evolution  system  was  integrated  using  the  Fehlberg  fourth- 
fifth  Runge-Kutta  method  based  routine  rkf45.f  over  the  time  interval  [0,30]  sec.  It 
should  be  noted  that  different  combinations  of  the  input  signal  u(t)  and  initial  conditions 
for  the  observer  state  were  also  tested  which  resulted  in  different  convergence  properties 
of  the  state  and  parameter  errors.  We  present  only  the  results  with  the  above  choice 
since  they  provided  the  required  level  of  persistence  of  excitation  needed  for  parameter 
convergence. 

The  evolution  of  the  state  components  x\{t)  and  X2 (t)  is  depicted  in  Figure  1  along 
with  their  estimates  over  the  interval  [—1, 10]  seconds.  Convergence  of  the  observer  state 
to  the  plant  state  can  easily  be  concluded.  The  evolution  of  the  parameter  estimates  along 
with  their  actual  values  are  depicted  in  Figure  2.  Parameter  convergence  is  also  achieved 
as  is  evident  from  the  figure. 

References 

[1]  R.  F.  Curtain,  Linear  operator  inequalities  for  stable  weakly  regular  linear  systems , 
manuscript,  (1997). 

[2]  R.  F.  Curtain,  M.  A.  Demetriou,  and  K.  Ito,  Adaptive  compensators  for  per¬ 
turbed  positive  real  infinite  dimensional  systems ,  manuscript,  (1997). 

[3]  R.  F.  Curtain  and  H.  J.  Zwart,  An  Introduction  to  Infinite  Dimensional  Linear 
Systems  Theory ,  Texts  in  Applied  Mathematics,  Vol.  21,  Springer- Verlag,  Berlin,  1995. 

[4]  K.  Ito  and  F.  K APPEL,  A  uniformly  differentiable  approximation  scheme  for  delay 
systems  using  splines ,  Applied  Mathematics  and  Optimization,  23  (1991),  pp.  217-262. 


0.5 


rii  and  rn(t)  xio_4  ri2andri2(t) 


Figure  2:  Evolution  of  the  estimates  r ij(t)  (solid)  and  their  actual  values  (dashed). 


[5]  H.  K.  Khalil,  Nonlinear  Systems ,  Macmillan,  New  York,  1992. 

[6]  K.  S.  Narendra  and  A.  M.  Annaswamy,  Stable  Adaptive  Systems ,  Prentice  Hall, 
Englewood  Cliffs,  NJ,  1989. 


